Effect of initial conditions on Glauber dynamics in complex networks 
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The effect of initial spin configurations on zero-temperature Glauber spin dynamics in complex 
networks is investigated. In a system in which the initial spins are defined by centrality measures 
at the vertices of a network, a variety of non-trivial diffusive behaviors arise, particularly in relation 
to functional relationships between the initial and final fractions of positive spins, some of which 
exhibit a critical point. Notably, the majority spin in the initial state is not always dominant in 
the final state, and the phenomena that occur as a result of the dynamics differ according on the 
initial condition, even for the same network. It is thus concluded that the initial condition of a 
complex network exerts an influence on spin dynamics that is equally as strong as that exerted by 
the network structure. 

PACS numbers: 89.75.Fb, 05.50.4-q, 64.60.Cn, 89.75.Hc 



The study of complex networks has generated broad in- 
terest, particularly in the field of nonlinear physics [U, 
Dynamic processes and the dependence of such processes 
on networks are an important topic in the field of complex 
networks, and many situations that result in non-trivial 
global phenomena have been discovered, despite the rela- 
tive simplicity of local interactions. The Glauber dynam- 
ics of the Ising model, involving a set of spin-like binary 
variables and corresponding local interaction, are some of 
the simplest dynamics that occur in such networks. The 
Ising model with Glauber dynamics has been investigated 
as a model of the ordering process in spin glasses and 
other systems Q. The interaction patterns of Glauber 
dynamics on regular d-dimensional lattices and complex 
network structures have been considered in a number of 
studies, and it has been shown that such structures pro- 
mote the development of characteristic ordering from a 
completely disordered state [1, H, [E S @ • 

Yet beyond the effect of structural complexity itself, 
complex networks also exhibit complexity related to the 
dynamics peculiar to complex structures. Previous of 
Glauber dynamics have predominantly focused on the 
ordering dynamics that emerge from the completely dis- 
ordered state, that is, the state in which the initial spins 
are randomly distributed with an initial fraction r of posi- 
tive spins. However, the vertices in complex networks are 
not interconnected uniformly, instead they form a het- 
erogeneous distribution of their degrees or locations. As 
a result of this heterogeneity, vertices in networks must 
be characterized topologically, such as by the measure 
of centrality Q, which represents a set of parameters for 
individual vertices that are determined according to net- 
work topologies. If vertices are discriminated by central- 
ity measures, an initial configuration can be considered 
that is not randomly distributed, but rather topologically 



biased. In such a case, it becomes possible to introduce 
some form of ordering into an initially disordered state. 
Given initial conditions, Glauber dynamics can also be 
considered as controlling the diffusion of spins from ar- 
bitrarily selected vertices in the initial condition. Such 
a semi-ordered initial condition, in addition to the effect 
of the network structure, is expected to have a signif- 
icant effect on the resulting dynamics. An attempt to 
demonstrate this effect through exploration of a simple 
two-state diffusion process has been reported previously 
[lol |. However, while the numerical results can be cate- 
gorized into several classes, they are considered merely 
preliminary, since the model is unfamiliar and numerical 
experiments alone are insufficient to confirm such cate- 
gorization. 

In the present study, the dependence of zero- 
temperature Glauber dynamics, one of the simplest forms 
of Glauber dynamics, on the initial spin configuration is 
investigated by Monte Carlo simulations of several types 
of complex networks. The dynamics are taken from the 
arbitrary, semi-ordered initial spin configuration deter- 
mined by the characteristics of the network structure. 

Zero-temperature Glauber dynamics are considered 
here for the case of a complex network with spin vari- 
ables (7 = ±1 located at the vertices of the network. The 
network is denoted by the adjacency matrix Aij , which 
is defined such that Aij = 1 if the vertices i and j are 
connected, and Aij = otherwise. In the present case of 
an undirected network, Aij is symmetric. The local field 
hi{T) acting on vertex i during time step r due to the 
spins of the neighboring vertices of the vertex i is given 

by 
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In the preceding study, similar dynamics were investi- 
gated as a two-state diffusion process [l^l by considering 
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the following model: 



(Tt{n + 1) 



sgn{hi{n)} if hi{n) ^ 
<j(n) if hi{n) — 0. 



(2) 



In this model, the spins of vertices are updated syn- 
chronously at each time step n, that is, the spins at all 
vertices are updated simultaneously as n progresses. In 
the present study, a general Monte Carlo method is cm- 
ployed for numerical simulations, and the spins of indi- 
vidual vertices are updated asynchronously. The rule for 
updating the spin of a randomly selected vertex is as fol- 
lows: 



sgn{/i,(T)} if K{t)^^ 

±1 (probability 1/2) if /ij(r) 0. 



(3) 

The progressive behavior of the spins is investigated using 
a normalized time step oit = t/N. 

Spins on rf-dimensional lattices have been investigated 
by a number of researchers, and it has been shown that 
despite the simplicity of zero-temperature Glauber dy- 
namics, non-trivial phenomena arise under such regimes, 
even for regular lattices [HE El- The Glauber dynam- 
ics of other types of complex networks, such as the Watts- 
Strogatz network [l4l| . have also been examined. Boyer 
and Miramontes [6], by analyzing the ordering dynam- 
ics of the Watts-Strogatz network, revealed a nonequi- 
librium ordering process induced by 'shortcuts' in the 
WS network, with occasional incomplete ordering. Ran- 
dom graph network of Erdos-Renyi and scale-free 
networks (i.e., Barabasi- Albert [l6|) have also been in- 
vestigated 0, for which Glauber dynamics do not lead 
to a fully ordered state. In such systems, the dynam- 
ics may become trapped in a set of partially ordered 
(meta)stable states, even when the system is finite. Zhou 
and Lipowsky described a mean-field approach for 
analyzing the dynamics in generic uncorrelated complex 
networks having arbitrary degree distributions. In scale- 
free networks with P{k) ~ a characteristic power- 
law exponent jc at which a dynamic transition takes 
place was identified. The voter model, which is similar to 
zero-temperature Glauber dynamics, has also been stud- 
ied in complex networks O, [H, HO, ISI . Remarkably, 
some of these studies have indicated that the majority 
population in the initial state may not dominate asymp- 
totically [21]. 

In these previous studies, the initial spins of the net- 
work vertices are considered to be randomly distributed, 
and it remains to be discussed whether a degree of order 
in the initial state affects the network dynamics. The 
measure of centrality has been used in network science to 
represent the characteristics of vertices relative to other 
vertices [1, [l^ . The use of centrality to characterize ver- 
tices thus allows a certain degree of order to be intro- 
duced into the initial spin configuration to which Glauber 
dynamics are applied. In the present study, the depen- 
dence of the dynamics on the initial spin configuration 



is examined using three kinds of centrality measures: de- 
gree centrality Cdeg{v), closeness centrality Cciose(v), and 
betweenness centrality Cbet{v) [II]- The clustering coef- 
ficient Cciust(v) is also considered as a characteristic of 



vertices 
as follows: 



141. These measures are defined for a vertex v 
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where fc.y is the degree of the vertex v, dcivji) is the 
distance between vertices v and t, ast and astiv) are the 
number of shortest paths from s to i and from s to t via 
V, and E is the number of edges between the neighbors 
of vertex v. 

Corresponding to the centralities of each vertex, the 
rN vertices with the largest centrality are assigned the 
positive spin state in the initial condition (cr(0) = -1-1), 
while the remaining {l~r)N vertices are assigned (t(0) — 

— 1. Here, r is the initial fraction of positive spins, and 
N is the number of vertices in the network. Four types of 
complex network models are investigated using these ini- 
tial conditions: Random graph network of Erdos-Renyi 
[m . and three types of structured scale-free networks 

- the Barabasi- Albert (BA) network [l^, the Klemm- 
Eguflz (KE) network (also reffered to as the highly clus- 
tered scale- free network) [l^, and the connecting near- 
est neighbor (CNN) network [2^. In the highly clus- 
tered scale-free network, the mixing parameter is set to 
/i = 0.1. All the scale- free network models are cate- 
gorized as evolving network models in which the distri- 
butions of the degree k follow a power-law distribution 
P{k) ^ k^'' . However, the networks differ in that the BA 
network has a low clustering coefficient C ~ N~^-''^ while 
the other two networks are much more highly clustered 
(C ~ C(l)), and the CNN network has a positive (as- 
sortativie) degree correlation. See References [l6l[2^[25| 
for further details. The same network models were used 
with a fixed number of vertices and average degrees in the 
preceding study [l^]- In the present paper, the scaling 
effects of system size are also considered. In the series of 
simulations described below, all the vertices in the net- 
works are confirmed to be connected by edges so taht they 
form a single component. 

The progression of positive spin fraction R+(t) at 
t — t/N, starting from fixed initial fractions of r = 0.4 
and r = 0.7 and a variety of initial orderings of positive 
spins is shown in Fig. [1] The results shown through- 
out this report are averages obtained over several simu- 
lations. In the random and BA networks, the dynamics 
reach a fully ordered state in which all of the vertices 
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FIG. 1: Temporal evolution of R+{t) in a random network (A), a BA network (B), a KE network (C) and a CNN network (D). 
Initial fractions r = 0.4 and 0.7. The different symbols denote the ordering of the initial state, according to degree centrality 
(□), closeness centrality (o), betweenness centrality (A ) and clustering coefficient (v). A'^ = 18000 and (k) = 10. 



share the same spin state. With r — 0.4 in the initial 
state, positive spins derived from vertices having larger 
centralities, despite being the minority, spread over the 
entire network in all cases except the BA network with 
the initial ordaring based on the clustering coefficient 
measure. In the KE and CNN networks, the dynam- 
ics reach a metastable state consisting of two coexisting 
spin states with fraction r', which differs from r. 

Figures [5] and [3] show the relationships between the 
initial fraction r = R+{0) and the convergence fraction 
r' = i?+(oo) of positive spins in systems having various 
sizes and average degrees. Several classes of character- 
istic functional forms can be observed. If the respective 
numbers of two spin states are not equal (i.e., r ^ 0.5), 
and the two states are randomly distributed, it is obvious 
that the state having the larger quantity in the initial 
state will become dominant at t = oo. This was veri- 
fied by performing preliminary simulations (Fig. [¥]). In 
the random and BA networks, the convergence fraction 
of the two spin states exhibits a step-like change with 
respect to initial conditions, where the critical fraction 
r = Tc < 0.5 is non-trivial. Notably, Vc is always lower in 



the BA network than in the random network. 

In contrast, the KE and CNN networks produce much 
more complex variation in the convergence fraction of the 
two spin states with respect to ordering of the initial spin 
configuration, reflecting the emergence of a metastable 
state. The initial states with ordering defined by degree 
and betweenness centrality exhibit convex or sigmoidal 
functions over r, whereas the relationship for initial or- 
dering defined by closeness centrality is almost linear. A 
critical point of r = rc can also be identified in these 
traces, around which r' jumps or rapidly increases in 
some cases. If the initial spins are completely disordered, 
Tc should take a value of 0.5. However, the value of Vc 
differs according to the initial condition applied. 

Figures [2] and [3] show the r-r' relationships for systems 
of various sizes and average degrees. Although some dif- 
ferences can be seen in the detailed shape of the rela- 
tionship and the value of comparison of Figs. [5] and 
[3] suggest that the class of the characteristic function is 
independent of the system size and the average degrees, 
except in the case of the CNN network. To investigate 
the scaling effects in more detail, the critical behavior of 
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FIG. 2: Convergence fraction of positive spins r' as a function of initial fraction r, derived from the initial ordering by (A) 
degree centrality, (B) closeness centrality, (C) betweenness centrality and (D) clustering coefficient, on the random network 
(□), the BA network (o), the KE network (A ) and the CNN network (V). Open symbols and closed symbols denote data for 
system size N = 18000 and iV = 36000, respectively, (fc) = 10. 




FIG. 3: Same quantities and symbols as Fig. O Open and closed symbols and filled symbols denote the plot for the average 
degree (k) — 10 and (fc) — 18, respectively. The system size is A'^ = 9000. 
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FIG. 4: Convergence fraction of positive spins r' as a function 
of initial fraction r derived from a completely disordered (i.e., 
randomly distributed) spin state. A'' = 18000, (k) = 10. 



the r-r' relationships was examined over a wider range 
of system size considering a representative set of initial 
orderings. A finite scaling on the CNN network can be 
recognized in both cases shown in Fig. \5\ As the sys- 
tem size increases, the critical point Tc asymptotically 
approaches a certain value, and the transition at Tc is 
not abrupt. No such scaling effects can be recognized 
on the KE network from Fig. [51 where the shape of the 
function and the value of Tc remains relatively constant, 
independent of system size. 

The characteristics of the metastable state emerging 
in the KE and CNN networks are illustrated using the 
fraction of vertices that flip spin (i.e., the vertices that 
connect to an equal number of different spins, or the ver- 
tices located at the boundary of a domain of identical 
spins), and the number of connected domains of positive 
spins ait = t/N. These parameters are denoted by ri^(t) 
and £>+ (t) , respectively. Figure [7] shows an example of 
the variation in nA(t) and D^{t) with r in systems with 
initial ordering defined by degree centrality. It can be 
seen that D+{t) does not always reach unity, suggesting 
that the metastable state can consist of several separated 
clusters of identical spins. Moreover, nA{t) does not con- 
verge to zero, indicating that the metastable state is not 
static but rather a stationary active state. The results 
are essentially the same for other initial orderings. It has 
been reported that a similar phenomenon might occur 
in a random network from a completely disordered ini- 
tial state [7| , implying that the final state consists of two 
large domains of opposite spin. A similar investigation 
was conducted using high-dimensional lattices The 
present simulations considering a range of complex net- 
work structures and arbitrary initial conditions suggest 
an alternative outcome is possible. 

The plots of nA{t) and D+{t) in the final metastable 
state (at t — oo) with respect to r for various initial 
conditions (Fig. [5]) indicate that the final states for the 
two networks are quite different. On the KE network, 
n^(cx)) exhibits a clear peak for each set of initial condi- 
tions, whereas a more complicated noncontinuous func- 
tion appears on the CNN network. For a completely dis- 
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FIG. 5: Convergence fraction of positive spins r' around the 
critical point as a function of initial fraction r for the CNN 
networks of various sizes with initial orderings defined by (A) 
degree centrality and (B) betweenness centrality. 



ordered initial configuration, nyi(oo) reaches a maximum 
at r = 0.5 and forms a symmetric function. For most 
initial conditions on both networks, the initial fraction 
r affording the maximum n^(oo) is consistent with the 
critical fraction r = Tc in Fig. [2l It is therefore safe to 
conclude that the metastable state arising from r ~ Tc is 
somewhat unstable. The number of connected domains 
of positive spins (Fig. [5]) reveals another characteristic. 
At small r, positive spins tend to shrink and become ex- 
tinct, thus £'+(oo) = 0. As r increases, several clusters 
of positive spin are able to persist, causing D^{oo) to 
increase. As r increases further, the fragmented domains 
begin to agglomerate to form fewer, larger clusters, until 
a singl large domain emerges (i.e., Z)+(oo) = 1). The de- 
gree of fragmentation depends on the initial conditions 
and the networks. It should be noted that D^{oo) is 
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always equal to or very close to unity in both networks 
when the order of the initial spin configuration is defined 
by the closeness centrality measure. It is suspected that 
this gives rise to the linear relationship between r and r', 
as shown in Fig. [DJB). 

In summary, the effect of the initial spin configuration 
on zero-temperature Glauber dynamics acting on various 



types of complex networks was examined. Through a se- 
ries of numerical analyses it was revealed that non-trivial 
diffusive behavior that depends on both the initial con- 
dition and the network structure occurs in such systems. 
In some cases, the final dynamics reached a metastable 
state involving two coexistent spin states with several 
connected domains of identical spins, the marginal ver- 
tices of which flip continuously. The functional relation- 
ship between the initial fraction r of positive spins and 
the final fraction r' can be categorized into one of sev- 
eral patterns, some of which exhibit critical point re- The 
patterns and criticality were found on networks of various 
system sizes and average degrees. A finite scaling effect 
with respect to system size was identified on the CNN 
network, whereas the other networks examined exhibited 
no such scaling effect. 

Although a clear understanding of the origin of the 
behavior revealed in this study has yet to be obtained, 
some observations can be made on the general character- 
istics of these phenomena. As the KE and CNN networks 
are known to be more highly clustered than the random 
and BA networks, the final metastable state may origi- 
nate from the topological cohesiveness of such networks. 
However, the characteristic functional relationships be- 
tween r and r', and the origins of the critical point rc, 
remain unclear. One possible feature that could be re- 
sponsible for these characteristics is the correlation in 
networks [1^. The CNN network is known to have a 
positive degree correlation [2^ , and hence if the vertices 
of higher degree are more likely to be interconnected, 
the positive spins originating from these vertices will not 
disappear, even for a small initial fraction. For other 
initial orderings, the correlation of centralities may also 
play a key role in the unique dynamics observed. The 
results of the present study may therefore lead to novel 
understandings of the underlying characteristics of in- 
dividual vertices measured thus far by centralities. It 
is also expected that further characteristics of network 
topologies remain to be revealed. The nonequilibrium 
spin dynamics exhibited by these systems do not appear 
to be comprehensively explainable based on topological 
scaling or conventional measures of the underlying com- 
plex networks alone. Nevertheless, the present study has 
revealed a strong dependence of spin dynamics on the 
initial conditions determined by the topological hetero- 
geneity of the network. 
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